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Future cosmological measurements should enable the sum of neutrino masses to be determined 
indirectly through their effects on the expansion rate of the Universe and the clustering of matter. 
We consider prospects for the gravitationally lensed Cosmic Microwave Background anisotropies and 
Baryon Acoustic Oscillations in the galaxy distribution, examining how the projected uncertainty 
of « 15 meV on the neutrino mass sum (a 4 (t detection of the minimal mass) might be reached 
over the next decade. The current la uncertainty of ~ 103 meV (PJanck-2015+BAO-15) will 
be improved by upcoming ‘Stage-3’ CMB experiments (S3-I-BAO-15: 44 meV), then upcoming 
BAO measurements (S3+DESI: 22 meV), and planned next-generation ‘Stage 4’ CMB experiments 
(S4-I-DESI: 15-19 meV, depending on angular range). An improved optical depth measurement 
is important: the projected neutrino mass uncertainty increases to 26 meV if S4 is limited to 
^ > 20 and combined with current large-scale polarization data. Looking beyond ACDM, including 
curvature uncertainty increases the forecast mass error by « 50% for S4-I-DESI, and more than 
doubles the error with a two-parameter dark energy equation of state. Complementary low-redshift 
probes including galaxy lensing will play a role in distinguishing between massive neutrinos and a 
departure from a w = — 1, flat geometry. 

I. INTRODUCTION 

A central goal in both cosmology and particle physics 
is to measure the mass of the neutrino particles. The 
neutrino sector is still poorly understood and the mecha¬ 
nism that gives rise to their mass is unknown. There are 
thought to be three active neutrino species, with mass 
differences measured through solar, atmospheric, reactor 
and accelerator neutrino oscillation experiments (for re¬ 
views see e.g., Gonzalez-Garcia and Nir [T], Maltoni et al. 

[2], Smirnov [3], Feldman et al. [4]). The results imply 
a minimum total mass of 60 meV in a normal hierarchy 
with two lighter and one heavier neutrino, or 100 meV in 
an inverted hierarchy with two massive neutrinos. 

Cosmology provides an indirect probe of massive neu¬ 
trinos [e.g., ISHH]. Massive neutrinos behave initially like 
non-interacting relativistic particles, and then later like 
cold dark matter. As such they affect the expansion rate 
of the Universe, compared to a pure radiation or pure 
matter component, as well as modifying the evolution 
of perturbations at early times. They also modify the 
growth of structure through a suppression of the cluster¬ 
ing of matter on scales that entered the cosmic horizon 
while the neutrinos were relativistic. 

The current indirect 95% upper limit from cosmo¬ 
logical data on the sum of the neutrino masses is 
TjUIi, < 230 meV from the Planck measurements of the 
Cosmic Microwave Background (CMB), combined with 
Baryon Acoustic Oscillation (BAO) measurements from 
the Baryon Oscillation Spectroscopic Survey (BOSS) 

[I^ [19] . The limit is Hrrii, < 680 meV from the CMB 
alone [H]. Tighter limits have been found including 
Lyman-a forest measurements from quasars in the BOSS 


survey (Etoj/ < 120 meV) [IHl HI], but the result de¬ 
pends on numerical hydrodynamical simulations which 
may contribute additional systematic uncertainty. 

Recent forecasts of mass limits for upcoming cosmo¬ 
logical datasets, including galaxy lensing and cluster¬ 
ing, redshift-space distortions, the kinematic Sunyaev- 
Zel’dovich effect, and counts of galaxy clusters, have 
been studied extensively (e.g., Font-Ribera et al. [22) . 
Villaescusa-Navarro et al. [53], Mueller et al. jH]), show¬ 
ing the promise of a wide range of future cosmological 
data to target a neutrino mass measurement. In this 
paper we focus on the combination of lensed CMB and 
BAO measurements, datasets which do not require de¬ 
tailed modelling of non-linear structure formation, or an 
understanding of galaxy bias. The gravitationally lensed 
CMB measures the growth of structure at times typi¬ 
cally before the Universe was half its current age, and 
on angular scales larger than « 100 Mpc, so is domi¬ 
nated by linear physics. Studies of this combination have 
been reported in Hall and Challinor [^, Abazajian et al. 
[53] , Wu et al. [57] and Pan and Knox [5S] , with a 4cr de¬ 
tection of neutrino mass forecast for the next generation 
of experiments. In this paper we investigate this further, 
exploring the dependence on experimental details and on 
parameter degeneracies. 

In [jn] we give a brief review of the cosmological ef¬ 
fects of neutrinos, and in §III| study how the mass mea¬ 
surement may be reached step-wise using data collected 
during the coming decade. In §IV| we investigate the 
dependence on experimental details, and in we ex¬ 
plore degeneracies with other cosmological parameters. 
We conclude in §VI[ 


2 


II. COSMOLOGICAL EFFECTS OF 
NEUTRINOS 

Standard Model neutrinos are initially relativistic, fol¬ 
lowing a thermal distribution after decoupling from the 
primordial plasma when the Universe had a temperature 
of around fcsT « 1 MeV. The neutrino temperature de¬ 
creases as the scale factor grows, until their rest-mass en¬ 
ergy dominates and they become indistinguishable from 
cold dark matter. For a neutrino of mass mv the tran¬ 
sition occurs at z « ^^^( eomev ) P^ - current limits 
indicate a transition epoch of 120 < z < 460 for a nor¬ 
mal mass hierarchy. 

This limit implies that the neutrinos were still rel¬ 
ativistic when the CMB decoupled, so they would be 
indistinguishable from massless neutrinos in the pri¬ 
mary anisotropies. However, higher mass neutrinos be¬ 
come non-relativistic sooner, which reduces the early- 
time Integrated Sachs-Wolfe (ISW) effect. This gravi¬ 
tational redshift of the CMB photons arises while the 
non-negligible radiation component causes the potentials 
of the density fluctuations to evolve mm and affects 
the anisotropies on scales around the first acoustic peak 

m- 

This effect is not sensitive to masses which remain 
relativistic until well after decoupling [TOl |30] , but fur¬ 
ther information comes from probes of later-time large 
scale structure measurements. Massive neutrinos inter¬ 
act weakly, allowing them to free-stream out of overden¬ 
sities while relativistic, so the growth rate of matter per¬ 
turbations inside the horizon is suppressed compared to 
a universe with only cold dark matter. For comoving 
wavenumbers k ^ fcpSj Hu et al. [3T] show that the sup¬ 
pression of the matter power spectrum today, P{k), is 
proportional to the sum of the neutrino masses: 

PumAk) - P^m^=o(.k) _ _Q Q- / 

■ 'vO.leVy 0.136 y ’ 

( 1 ) 

where the comoving free-streaming scale is given by 

as illustrated in Fig. [^and e.g., [26], for models with fixed 
total matter density. For current limits this scale is esti¬ 
mated to lie in the range 0.005 < kps ^ 0.011 [h Mpc“^]. 

The suppression of small-scale power can be probed 
using galaxy clustering and the gravitational lensing of 
galaxies. These are promising avenues for neutrino mass 
measurements [e.g., [55], although these observables are 
sensitive to non-linearities in the matter power spectrum 
and scale-dependent galaxy and shape biases [35] . An al¬ 
ternative route is through the gravitational lensing of the 
CMB (see e.g., [33] for a review). Here the CMB photons 
are deflected by the large-scale structure, integrated over 
the photon path since decoupling. 

Following [33], the CMB convergence angular power 
spectrum, Cf”, is a weighted projection of the matter 


power spectrum P{k,x)\ under the Limber approxima¬ 
tion, 



where xh is the comoving horizon size, fk (x) relates line- 
of-sight comoving distances and transverse comoving dis¬ 
tances in a curved universe, and the window function 
W (x) is 

for X < X* s-ud zero otherwise. Here a(x) is the scale 
factor and x* is the radial comoving distance to the last¬ 
scattering surface. This power spectrum is sensitive to 
as shown in Fig. and does not depend on galaxy 
bias or detailed non-linear modeling. In practice it is 
reconstructed from CMB temperature and polarization 
maps using a four-point function [e.g.. 134] . 

The CMB temperature and polarization angular power 
spectra, {Cf^,Cj ^are also modified by 
lensing, which smears the acoustic peaks by adding vari¬ 
ance to the apparent scale of a mode, converts A-mode 
polarization into H-mode polarization, and adds small- 
scale power in T, E and B [e.g.,|33]. The approximate 
effect of massive neutrinos is shown in Fig. [^for the E- 
mode polarization, where we artificially amplify the ef¬ 
fects of neutrinos on the CMB lensing, rather than the 
primary CMB, by varying the amplitude of the lensing 
potential. Increasing the neutrino mass has a similar ef¬ 
fect to decreasing the lensing amplitude. Compared to 
the power spectra, the reconstructed convergence field 
contains more information on the neutrino mass [35] . 

Massive neutrinos also affect angular diameter dis¬ 
tances dA{,z) and the expansion rate H(z), as their evolu¬ 
tion differs from a pure radiation or pure matter compo¬ 
nent [e.g., 155 ]. These can be measured using a ‘standard 
ruler’ method that is relatively free of systematic uncer¬ 
tainties: the primordial oscillations in the photon-baryon 
fluid are imprinted in the galaxy distribution as Baryon 
Acoustic Oscillations (BAO). The comoving scale of the 
oscillations is fixed by the sound horizon at decoupling, 
Tg, which is not significantly affected by neutrino masses 
given current limits and is in the linear regime of density 
perturbations (« 150 Mpc). The observed spherically- 
averaged BAO angular scale for galaxies at redshift z is 
sensitive to the parameter combination rg/dviz); dy is 
the volume distance [ 55 ] . 

dviz) = [cz(l-I-z)^(i^(z)iL“^(z)] . (5) 

For fixed cold dark matter density, more massive neu¬ 
trinos increase the total late-time non-relativistic matter 
content, which increases the volume distance, as shown 
in Fig [2 
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FIG. 1: Effect of neutrino mass on CMB power spectra and BAO distance scales. Top: Fractional change of the matter power 
spectrum today P{k) (left) and CMB convergence power spectrum C'f'' (right) with neutrino mass Tirriv, for fixed physical 
dark matter density, + Tlvh^■ Suppression of power is due to neutrino free-streaming. Bottom left: Lensed CMB E-mode 
power spectrum with varying amplitudes of the lensing potential Aiens, approximating and exaggerating the effect that massive 
neutrinos have on the CMB polarization spectrum. Bottom right: BAO distance ratio rsjdy for fixed Qa and Massive 

neutrinos behave like additional matter in the BAO redshift range, decreasing Ho and increasing the volume distance dy- 


III. IMPROVEMENTS IN THE NEXT DECADE 

We consider how upcoming and planned CMB and 
BAO experiments will improve the current limits on the 
sum of the neutrino masses, building on previous analy¬ 
ses [SaEMSH]. 

A. New data 

The current state-of-the-art for the CMB is the Planck 
and WMAP satellite data, including the first analysis of 
the full-mission Planck data [El [37]. Improved small- 
scale CMB measurements are currently being made by 
the ‘Stage 2’ ground-based experiments: ACTPol, SPT- 
Pol and POLARBEAR [38b40j . These will soon be up¬ 
graded to ‘Stage 3’ (hereafter, S3) with new detectors 
and sensitivity to multiple frequencies. Here we con¬ 
sider an S3 ‘wide’ experiment that maps 40% of the sky. 


and a ‘deep’ experiment that maps 6% of the sky. The 
‘wide’ experiment is similar to AdvACT specifications 
[41] and ‘deep’ to SPT-3G [42]. These experiments are 
expected to take data during 2016-19, and the specifica¬ 
tions we adopt are given in Table We also anticipate 
data from S3 experiments targeting larger angular scales 
(e.g., CLASS), but do not consider these specifically. In 
addition, we can expect a complete analysis of the Planck 
polarization data, optimistically including reliable large- 
scale polarization data. 

Beyond S3, a ‘Stage-4’ (S4) experiment - or set of ex¬ 
periments - is being developed by the CMB community 
that may cover at least half the sky to typical noise levels 
of 1 /rK-arcmin |E] • There are also proposed space-based 
experiments including LiteBIRD and PIXIE |43l 1^ . 
and we approximate their role with a cosmic variance- 
limited large-scale polarization measurement (‘CV-low’) 
that could supplement S3 ground-based data. 

On the BAO front, current state-of-the-art measure- 
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FIG. 2: Top and middle: Fractional change in the conver¬ 
gence ft (top) and E-mode (middle) power spectrum with 
neutrino mass, for fixed with expected uncer¬ 

tainties for S3-wide and S4 CMB data. A higher neutrino 
mass has less lensing, decreasing the E-mode peak smooth¬ 
ing. Bottom: Fractional change in distance ratio rsjdv, with 
uncertainties from current (BAO-15 [18]) and forecast (DESI, 
|22|i BAO data. Here is fixed. 


Experiment /sky 

Beam 

(arcmin) 

AT 

(pK-arcmin) 

AP 

(/rK-arcmin) 

S3-wide 

0.4 

1.4 

8.0 

11.3 

S3-deep 

0.06 

1 

4.0 

5.7 

S4 

0.4 

3 

1 

1.4 

CV-low 

0.4 

60 

1 

1.4 


TABLE I: Upcoming (S3) and proposed (S4, CV-low) CMB 
experiments that will add to Planck. The Planck specifica¬ 
tions we use are in the Appendix. 


ments come from the BOSS ‘LowZ’ and ‘C-MASS’ galaxy 
samples at z = 0.32 and z = 0.57 [18]. These are supple¬ 
mented by data from the Six Degree Field Galaxy Red- 
shift Survey at z = 0.11 [IS] and the SDSS MGS sample 
at z = 0.15 [IS]. Improved measurements are being made 
by the eBOSS survey which will survey a deeper sample 
m- A significant advance should be made with the DESI 
spectroscopic survey, due to begin in 2018, which is ex¬ 
pected to measure the BAO distance ratio from redshifts 
0.15 < z < 1.85 in bins of width Az = 0.1 to percent- 
level precision [HI SH] • 


B. Forecasting methods 


We use a Fisher-matrix forecasting method to predict 
the neutrino mass uncertainties. For a model defined by 
parameters 6 the expected Fisher matrix is 


MS) 


\np{9\d) \ 

dOidOj / ’ 


( 6 ) 


where p(0|d) is the posterior distribution for 0 given data 
d. The forecast parameter covariance is then given by 
the inverse of the Fisher matrix, Here our 

data are the lensed TT, TE, and EE GMB power spectra, 
reconstructed GMB convergence power spectrum kk, and 
BAO distance ratio measurements rg/dy{z). Our param¬ 
eters are the standard six ACDM parameters, plus the 
neutrino mass sum, as well as possible extension param¬ 
eters including curvature and dark energy. Our methods 
are summarized in the Appendix, including choices made 
about the fiducial model, choice of parameter basis, and 
step-size for calculating derivatives. We also describe val¬ 
idation of our numerical code. 

We use the lensed CMB power spectra and the conver¬ 
gence power spectrum as our CMB observables, which 
differs from the approach in [HI , but more closely 

follows the ‘real’ data analysis: the CMB sky we see is 
lensed, and it is a difficult inverse problem to infer the 
unlensed sky [e.g.,|35]. Using unlensed spectra in fore¬ 
casts removes information contained in the lensed tem¬ 
perature and polarization fields. However, it is challeng¬ 
ing to construct the full covariance matrix for the lensed 
power spectra and convergence power spectrum: the T, 
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Q and U fields are all lensed by the same leasing poten¬ 
tial, which correlates the power spectra and adds addi¬ 
tional non-Gaussian covariance. This is explored in de¬ 
tail in [Ml HO]- In this analysis we make the approxima¬ 
tion of discarding BB information, and assuming Gaus¬ 
sian uncorrelated errors in TT, TE, EE, and kk. This 
is likely a good approximation for S3 data, but could 
under-estimate certain parameter errors for S4-type data 
by up to « 20% |M]- We use CAMB for evaluation of all 
relevant CMB power spectra [51j . 

For the noise levels of Planck, we consider two cases, 
^Planck-2015’ (P15) that produces cosmological con¬ 
straints which closely match the published results m, 
and ‘Planck-poV which includes TE and EE data com¬ 
ing from the polarization measurements of the High- 
Frequency Instrument (HFI), including large-scales. The 
specifications are given in the Appendix. For Planck-pol, 
noise levels are approximated by taking temperature sen¬ 
sitivities from P15 and assuming the per-channel noise 
scalings from temperature to polarization in the Planck 
Blue Book |M] • This is likely to be over-optimistic at the 
largest scales. 


For the GMB power spectra, we set a maximum mul¬ 
tipole for the recoverable information: = 3000, 

^max = 4000 and = 3000 for the future S3 and S4 ex¬ 
periments. Smaller scales are likely hard to extract due 
to extragalactic foreground contamination. We assume 
white noise, and do not include additional foreground 
uncertainty beyond the multipole cuts outlined above, 
although the expected S3 white noise level includes some 
foreground inflation |41] . We also set a minimum multi¬ 
pole of fmin = 50 for S3 due to the challenge of recovering 
large-scales from the ground, and consider two options for 
S4: £min = 50 and = 5. For S3 and S4 we include 
Planck data for 2 < i < ^min, and our nominal analysis 
uses ‘Planck-poV unless stated otherwise. We consider 
the importance of the large-scale polarization measure¬ 
ments in §IV A| 


We use the quadratic-estimator formalism of Hu and 
Okamoto |S3] to calculate the GMB convergence noise 
spectrum N^'^. This uses the coupling of otherwise un¬ 
correlated modes in temperature and polarization to re¬ 
construct the lensing potential. Iterative delensing pro¬ 
cedures are able to reduce the effective noise level of 
the estimated lensing field, particularly for the low-noise 
(AP < 4^K-arcmin) future experiments considered here 
[54] . We consider the impact of this process in 


illlC 


For the BAG measurements we use the published un¬ 
certainties on the distance ratio Tsjdy for the current 
BAG data described in ? HI [TSl[Ml@5] , labelled as ‘BAG- 
15’. For DESI we use the forecast uncertainties on dA{z) 
and H{z) given in Font-Ribera et al. [22] to estimate 
the expected r^ldy uncertainties, summarized in the Ap¬ 
pendix. We do not use broadband shape information in 
the galaxy power spectra. 



FIG. 3: Forecast marginal posterior constraints on the sum of 
the neutrino masses Em,/ within a ACDM-|-Sm„ model, as¬ 
suming Gaussian error distributions. The current uncertain¬ 
ties (P15-I-BAO-15) are expected to improve rapidly, with S3 
GMB data and DESI BAG data expected by ~2020. 



Em,/ [meV] 


FIG. 4: Expected joint constraint (68% GL) on the neutrino 
mass sum Em,/ and physical cold dark matter density Vljp 
within a ACDM-I-Em,/ model. The BAG constraint, sensi¬ 
tive to the total late-time cold dark matter density, is almost 
orthogonal to the GMB lensing constraint, breaking the de¬ 
generacy. 
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C. Expected constraints 


We find forecast marginal Itr uncertainties on the sum 
of the neutrino masses T,mi, of 


(T(Emy) 

meV 


'l03 (P15 + BAO-15) 

44 (S3-wide + BAO-15) 

< 22 (S3-wide + DESI) 

19 (S4 (e > 50) + DESI) 

15 (S4 (^ > 5) + DESI) 


( 7 ) 


where S3 and S4 include Planck-pol at large-scales. If 
we replace S3-wide with the S3-deep survey we find 
= 53 meV combined with BAO-15, and 25 meV 
with DESI. The expected constraints are summarized in 
Fig. i and are consistent with findings in Wu et al. m- 
These forecasts imply that if the neutrino hierarchy is 
inverted, with mass sum > 100 meV, we may have > 2a 
evidence for non-zero mass in the next few years from S3 
data, and an almost 5cr detection in years with DESI. 
If instead it is the minimal mass of 60 meV, a 2 — 3cr indi¬ 
rect measurement should be reachable in five years, with 
stronger evidence from the subsequent experiments. 

As illustrated in Fig. there is a strong positive cor¬ 
relation between the neutrino mass and cold dark matter 
density in the CMB observables. This arises predomi¬ 
nantly from the competing influence of these parameters 
on the lensing signal. Increased neutrino mass suppresses 
small-scale power, while increasing cold-dark matter den¬ 
sity boosts small-scale power by shifting matter-radiation 
equality to earlier times; this shortens the radiation- 
domination epoch in which sub-horizon modes of the 
gravitational potential decay, enhancing the small-scale 
amplitude of structure [5SJ [SS]. Conversely, BAO 
data constrain the sum of the CDM and massive neu¬ 
trino density, since it is this combination that affects the 
angular diameter distance and expansion rate. 

The power of the BAO and CMB data lie in their com¬ 
bination. BAO measurements alone cannot constrain the 
primordial parameters As and Ug, nor the optical depth 
to reionisation r; the constraints lie in the 4-dimensional 
subspace spanned by and 9a (being 

most constraining in the directions corresponding ap¬ 
proximately to and 9a)- S4-alone could achieve 

a neutrino mass error of aCEm^) = 53 meV; with DESI 
BAO this is expected to tighten to a{Yjm^) = 19 meV, 
enough for a 3cr detection of the minimal mass. 

Comparing to previous results, the S4(£ > 5)-|-DESI 
forecast matches previous results and the ‘cur¬ 

rent data’ P15 -I- BAO-15 forecast is compatible with 
the published result: a{T,m^) < 0.23 eV at 95% confi¬ 
dence [ 15 ], with the forecast errors on other ACDM pa¬ 
rameters also matching closely. Here we chose a fiducial 
neutrino mass sum of 60 meV. We find a smaller uncer¬ 
tainty for higher neutrino mass (as in [25] 1. but the ef¬ 
fect is small when BAO data are included. For S4 alone, 
we find a{Yjmy) = 53 (41) meV for a fiducial mass of 


60 (120) meV, but with DESI this difference reduces to 
a{T,m^) = 19 (18) meV. We describe further tests of the 
code in the Appendix. 

We consider the impact of using an improved, 
likelihood-based lensing estimator (going beyond the 
first-order quadratic estimator of Hu and Okamoto [55] 1 
to reduce the effective noise power A^/^” of the lensing re¬ 
construction. Hirata and Seljak m show that Nf''^ can 
be reduced by a factor of two for an S4-like experiment; 
under this modification, we find only a 3% tightening of 
the neutrino mass constraint for S4-I-DESI. The improve¬ 
ment is small due to the contribution of cosmic variance 
(CV) to the lensing power spectrum uncertainties, and 
the degeneracy of neutrino mass with other ACDM pa¬ 
rameters; the S4 lensing power spectrum derived from the 
quadratic estimator is already CV-limited out to £ « 800. 

As in Benoit-Levy et al. [3S|, we find that there is useful 
information in the lensed power spectra. We show the 
relative impact on forecasted constraints for S3-wide and 
S4 in Table |H] (including Planck-pol at large scales but 
no BAO). Lensing reduces the neutrino mass uncertainty 
by a factor of « 6 compared to the unlensed CMB. As 
such, the two-point functions will provide increasingly 
more information as E-mode polarization measurements 
improve. Checking for consistency between the lensed 
observables will be an important systematic test for new 
CMB data. For example, marginalizing over the lensing 
effect in the 2-point functions, by introducing a variable 
lensing amplitude parameter Aiens jSS] , would isolate the 
impact of neutrinos on the 4-point function. 



Unlensed Lensed TT, TE, EE 

Unlensedd-KK 



(2-point only) 

(4-point only) 

S3, (j(Em,^): 

435 

75 

61 

S4, (j(Em^): 

363 

64 

53 


TABLE II: Impact of lensing on the neutrino mass constraint 
(in units of meV). Constraining power comes from both the 
lensed spectra (2-point) and the reconstructed lensing poten¬ 
tial (4-point). Gaussian uncorrelated errors are assumed. 


IV. DEPENDENCE ON EXPERIMENTAL 
DETAILS 

Since future CMB experiments are currently under de¬ 
velopment, we investigate the importance of certain ex¬ 
perimental details on the mass constraint. 


A. Importance of the reionization bump 

The amplitude of primordial power, Ag, is partially 
degenerate with Stoj/, since Ag increases the amplitude 
of clustering at small-scales, and decreases it. The 
amplitude Ag is not well determined by the primordial 
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CMB temperature anisotropy; an increased optical depth 
to reionization lowers the signal such that the normaliza¬ 
tion of the anisotropy measures the combination Ase~^'^ 
[55] , This leads to a degeneracy between r and Srui, 
which can be broken with precision measurements of the 
reionisation bump at multipoles ^ < 20 in polarization 

[nisn]- 

Here we explore the importance of making a robust 
optical depth measurement, considering three cases for 
S4: current WMAP measurements EH EB, optimistic 
future Planck-pol measurements (see Appendix), and a 
future S4 measurement that reaches the largest scales 
(^min = 5). We find forecast constraints of 


cr(STOi,) 

meV 


(S4 (£ > 50) 
(S4 (£ > 50) 
(S4 (£ > 5) 


WMAP-pol 

Planck-pol 

DESI) 


DESI) 
DESI) 

( 8 ) 

with the uncertainty on t reducing from 0.008 to 0.005 
to 0.003, respectively. This is compared to 0.013 for 
WMAP-pol from EE alone, i.e., improved CMB leasing 
data helps constrain t even when the neutrino mass is 
varied. Fig. shows the expected correlation between r 
and neutrino mass. 

Fig. [5] also shows the impact of reducing the mini¬ 
mum multipole of the S4 experiment on the neutrino 
mass constraint, supplemented with Planck-pol or the 
current WMAP-pol at the largest scales. There is a lim¬ 
iting plateau for S4 at fmin > 20, and a clear improve¬ 
ment as the polarization is better measured at increas¬ 
ingly large scales. The S4(£ > 5) -I- DESI limit reaches 
the cosmic-variance (CV) limit for CMB data^. 

We then also consider the relative importance of mak¬ 
ing a higher sensitivity small-scale measurement, versus a 
new large-scale polarization measurement. We start with 
an S3-type £> 50 experiment, and then either increase 
the ^ > 50 sensitivity, or supplement it with a new ‘CV- 
low’ large-scale measurement at ^ < 50. We find forecast 
constraints of 


meV 


'22 (S3 (£> 50) + Planck-pol-t DESI) 

= { 19 (S4 {£ > 50) -t Planck-pol + DESI) 

17 (S3 {£ > 50) -t CV-low -t DESI) 

(9) 

This indicates that a cosmic-variance-limited measure¬ 
ment of optical depth could be more valuable than more 
sensitive small-scale data, especially given that Planck- 
pol large-scale polarization data is itself not yet demon¬ 
strated to be free of systematic errors. 

We note that 21cm experiments, which map the bright¬ 
ness temperature of neutral hydrogen as a function of 


^ Pan and Knox m found that an S4 experiment combined with 
a CV-limited BAO experiment could tighten the neutrino mass 
constraint further, to cr{Ylmi,) = 11 meV. 



0 10 20 30 40 50 


^min 

FIG. 5: Top: The neutrino mass is correlated with the 
optical depth to reionization r (forecast 68% CL). Current 
data at ^ < 20 (WMAP-pol) would leave a degeneracy be¬ 
tween Eruiy and r that could be broken with improved large- 
scale polarization data. Bottom: the expected nentrino mass 
constraint as a function of the minimum multipole accessible 
to S4, indicating the benefit of reaching large scales. 


redshift, will probe the epoch of reionization |62j : combi¬ 
nation of this information with CMB-t-BAO would break 
the degeneracy and improve the neutrino mass 

constraint (investigated in Liu et al. [63)1. 


B. Importance of sensitivity and angular range 

For the particular case of an £ > 50 experiment cover¬ 
ing 40% of the sky at 3 arcmin resolution, combined with 
Planck-pol, we vary the white-noise sensitivity. The fore¬ 
cast neutrino mass limits are shown in Fig. |^for CMB- 
only, CMB-bBAO-15, and CMB-fBAO-DESI. There is 
clearly an improvement as the noise is reduced, and a 
significant gain is expected over current Planck measure¬ 
ments, but below white noise levels of « 10 /xK-arcmin 
there does not appear to be a substantial gain (as also 
seen in Wu et al. m)- 
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FIG. 6: The dependence of the nentrino mass constraint 
a(Em^) on CMB map sensitivity, for a 3 arcmin resolution 
experiment covering 40% of sky. 


It is not yet certain whether this « 10 /rK-arcmin noise 
level, over half the sky, will be achieved in practice from 
the upcoming S3 CMB experiments, or whether the lens- 
ing reconstruction will achieve the expected noise levels. 
Atmospheric, ground and foreground emission are typi¬ 
cal contaminants that would increase the effective noise 
in the maps and in the lensing reconstruction. New data 
from the current S2 experiments will help clarify the im¬ 
pact of non-white-noise on the lensing noise performance. 

Here we have continued to restrict our analysis to 
‘clean’ scales at £ < 3000 in TT, TE and £ < 4000 in 
EE. A foreground-free S3 or S4 experiment would contain 
information from I > 3000; the TT spectrum is signal- 
dominated to small scales {£ « 4800 for S4) where there 
are a large number of modes. However, uncertainty about 
extragalactic foregrounds will likely make this informa¬ 
tion inaccessible. We find that including £ > 3000 scales 
would tighten the neutrino mass error by « 10%. 

We restrict the modes available for reconstruction of 
the lensing potential to these same scales, and to £ > 50, 
since reconstruction on the largest scales has yet to be 
demonstrated (there are difficulties such as mean-field 
subtraction for masked fields [e.g., EH)- We find that 
the information, quantified by the term in the 

Fisher matrix, is concentrated in the multipole range 
100 < £ < 1000 for S3 and 200 < £ < 2000 for S4. This 
also means that detailed non-linear modeling should not 
be required, since at redshift z « 2 at the peak of the 
CMB lensing kernel (Eq. , the information peak corre¬ 
sponds to physical scales > 200 Mpc. 


V. HOW UNIQUE IS THE MASSIVE 
NEUTRINO SIGNAL? 

The measurable effect of neutrinos can be partly mim¬ 
icked by changes in other cosmological parameters, as ex¬ 
plored in e.g., Font-Ribera et al. [52], Benoit-Levy et al. 


[35] . Hamann et al. [5S]. Massive neutrinos affect the ex¬ 
pansion rate and angular diameter distance, but changes 
in curvature, dark energy history or the Hubble parame¬ 
ter can compensate to keep the well-constrained acoustic 
peak positions and structure essentially unchanged. 

Here we consider changes in the spatial curvature, and 
changes in the energy density of dark energy with time. 
Dark energy is invoked to explain the acceleration of 
the universe, and may take the form of a cosmologi¬ 
cal constant with constant energy per unit proper vol¬ 
ume, although an evolving dark energy equation of state 
{w yf —1) is not ruled out by current observations [19] . 
We consider two parameterizations for dark energy: these 
are the usual Taylor expansion in the scale factor for dy¬ 
namical dark energy [66] : 

w{a) =Wo + Wa{l-a), (10) 


with two free parameters Wg and Wa, and also the Doran 
and Robbers m model for early dark energy. 


flA{a) 


Da - He (1 - 

Ha -f 


+ He (1-0-3-°), (11) 


with parameters He and wg. This has a background ex¬ 
pansion similar to a massive neutrino for periods of the 
evolution of the universe. Previous work has considered 
neutrino mass constraints within this model [e.g. . 1681169] . 
We use the CAMB PPF module m and a modified ver¬ 
sion of CAMB from Calabrese et al. [68] to compute the 
power spectrum within these models. 

We take as our baseline the S4 (£ > 50)-|-DESI ex¬ 
periment. Marginalizing over simple extensions to the 
ACDM-|-Em^ model, we find 


meV 


'19 (ACDM-kSm^,) 

30 (ACDM-kSm^-kHfc) 

27 (ACDM-|-Stoi/-I-wo) 

46 (ACDM-l-STOi^-l-wo-l-Wa) 

37 (ACDM-kSTOiz-l-He-l-wo) 

64 (ACDM-kSTOiz-l-wo-l-rCa+Hfe) 


We discuss these parameter degeneracies in the remain¬ 
der of this section. We find that marginalizing over neu¬ 
trino number (and other extension parameters that 
modify the primordial CMB spectrum such as a running 
spectral index) have a < 10% effect on the predicted neu¬ 
trino mass uncertainties. These findings are summarized 
in Fig. [^ which also includes the corresponding S3 fore¬ 
casts. 

At first sight, the degradation of the neutrino mass 
estimate in the case of varying wg, Wa and H^ simulta¬ 
neously appears severe, more than tripling the error bar. 
However, this model has three extra parameters com¬ 
pared to ACDM, and within the Bayesian framework, 
one can rigorously ask whether additional parameters are 
required by the data, quantifying the trade-off between 
improving the fit against an increased complexity of the 










9 


80 


W 


S3(wide)+BOSS 
S3(wide)+DESI 
S4(£ > 50)+DESI 
S4(£ > 5)+DESI 



+ 


FIG. 7: Neutrino mass constraints forecast for different data combinations and simple one- or two-parameter extensions to the 
ACDM-fEm^ model (all except S4(.^ > 5) include Planck-pol information at low-^). With CMB-I-BAO data, neutrino mass is 
not correlated with the number of neutrino species, Neff, but is partly correlated with spatial curvature, fik, and with the dark 
energy equation of state, wo- Expected confidence levels are shown assuming the minimal total neutrino mass Em„= 60 meV. 


model [e.g., EH for discussion in the context of cosmolog¬ 
ical data]. This model selection approach would quantify 
the need for additional extension parameters, and would 
disfavor an over-parameterized model if it is not required 
by the data. In practice our challenge is likely to lie 
in distinguishing between different one-parameter exten¬ 
sions to ACDM: are we seeing non-zero neutrino mass, or 
could it be similarly well-explained by a small amount of 
curvature, or a small deviation from a w = —1 equation 
of state? 


A. Physical degeneracies 

To help understand these degeneracies, we note the 
Hubble parameter H{z) is given by 

2 

+ z)'^ + r2A(2'), 

at times after the neutrinos become non-relativistic, 
where Ha( 2 ) is the dark energy density, constant for 
w = —1, fir is the radiation density today (e.g., pho¬ 
tons and massless neutrinos), flrn = flc + + flu is the 

matter density today (CDM, baryons and massive neu¬ 
trinos) and ft^h"^ = / 93 eV is the physical massive 


H{z) 

Ho 


neutrino density today. The angular diameter distance 
is given by 

f 73?!^ sin ^yZ—ftk /q h{z') ) < 0, 

= ^.= 0 , 

[^sinh(A/n^/„^f^) Hfc>0. 

(13) 

Considering curvature, photons propagating in a non¬ 
flat universe follow curved geodesics, changing the angu¬ 
lar diameter distance to an object of fixed proper size, 
at a given comoving distance, relative to a flat universe. 
Varying curvature shifts the angular scale of the acous¬ 
tic peaks, so to remain consistent with the CMB data, 
the matter density and Hubble constant Hq must vary 
to keep the peak structure unchanged. This is the well- 
known geometric degeneracy m- 

This degeneracy is partially broken by CMB lensing 
measurements [e.g., E3], which are sensitive to the growth 
of structure in the late-time universe and therefore to the 
matter density flm and dark energy density Ha. For a 
fixed CMB acoustic peak scale, the effect of decreasing 
the curvature parameter H^, moving to a more closed 
universe, is to decrease the Hubble constant and increase 
flm- This enhances the amplitude of the CMB lensing 
power spectrum, as illustrated in Fig. which can be 
compensated by increasing Em,y. This leads to an anti- 
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FIG. 8: CMB convergence power spectrum for varying 
with other parameters holding the primary CMB fixed. De¬ 
creasing Qk requires a smaller Hubble constant Hq and in¬ 
creased growth rate. This has a similar effect to decreasing 
the neutrino mass. 


correlation between and 17^ when using CMB mea¬ 
surements alone, as shown in Fig.[^ 

The BAO constraint in the plane is approxi¬ 

mately orthogonal to the CMB-only constraint, because 
decreasing ftk decreases the volume distance to a given 
redshift (Eq. . This can be compensated by a smaller 
matter density, lowering the neutrino mass. These data 
are powerful in combination, and for S3-I-BOSS the neu¬ 
trino mass constraint is independent of curvature. How¬ 
ever, the neutrino mass constraint from S3-I-DESI or 
S4-I-DESI is expected to degrade by « 50% when allow¬ 
ing for curvature, as illustrated in Fig. 

Similar arguments apply to dark energy, which mod- 
ihes the background evolution according to its equation 
of state, but does not contribute to clustering. Increas¬ 
ing the dark energy equation of state parameter wq leads 
to an increased expansion rate, shifting the angular scale 
of the acoustic peaks, so to remain consistent with the 
CMB data the Hubble constant Hq decreases to keep the 
peak structure unchanged. Similar to the curved model, 
this increases the clustering, which can be compensated 
with larger neutrino masses. This gives the positive cor¬ 
relation between wq and T,m,j in CMB data, illustrated 
in Fig. [^and reported in e.g., ISSIITI]. The BAO degener¬ 
acy is also positively correlated however, so the neutrino 
mass uncertainty is inflated more than when allowing for 
curvature. Increasing increases the contribution of 
neutrinos to 0^, requiring a smaller Ha (in a flat uni¬ 
verse); the volume distance to a given redshift, and hence 
the BAO peak position, can then be preserved by increas¬ 
ing Wq. The early dark energy density parameter He, is 
anti-correlated with neutrino mass, due to their similar 




FIG. 9: Top; Forecast joint constraint on the neutrino mass, 
Em,/, and spatial curvature, Hfc, within a ACDM-|-Em,/-|-Hfe 
model. The BAO data breaks the anti-correlated degeneracy 
in the CMB data. Bottom: Forecast constraint on Em,/ and 
the dark energy equation of state, wo, marginalized over the 
ACDM and Wa parameters. 


effects on the background expansion [SHI IM] • ^ 

We find that the neutrino mass constraints from 
S3-I-DESI or S4-I-DESI are degraded by more than a fac¬ 
tor of two when allowing for a non-minimal (w yf — 1) 
dark energy equation of state in the cosmological model, 
as shown in Fig. Distinguishing a non-zero neutrino 
mass parameter from a universe with dark energy be¬ 
yond the cosmological constant will be difficult with the 
CMB-I-BAO data combinations considered here. Allow¬ 
ing for freedom in both dark energy and curvature (i.e. 
a ACDM-|-Em,/-|-Hfc-|-wo-l-'iUa model) degrades the con¬ 
straint further to (jiYim^) = 64 meV for S4-1-DESI, but 
would include three new parameters. 


^ Here we adopted a fiducial He = 0.007, allowed by current data 
CD- A smaller Qe would improve the neutrino mass constraint, 
as the parameters are anti-correlated and Qe cannot be negative. 
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FIG. 10: Galaxy leasing power spectra for two models degen¬ 
erate in GMB spectra and BAG distance ratios. They could 
be distinguished using a future galaxy weak leasing survey 
(/sky = 0.5, rigai = 10 arcmin~^, and source redshift distribu¬ 
tion n{z) as for the LSST ‘gold sample’ [76]). 


B. Breaking degeneracies with complementary 
measurements 

We have focused so far on the minimal combination 
of future CMB and BAO data. To break these re¬ 
maining degeneracies between neutrino mass and curva¬ 
ture/dark energy parameters, we would turn to other 
large-scale structure probes such as the galaxy power 
spectrum, redshift-space distortions, galaxy weak leas¬ 
ing, the kinematic Sunyaev-Zel’dovich effect and galaxy 
cluster counts, which measure the growth of structure at 
later times, and Type la supernovae which better con¬ 
strain the expansion rate. 

We have not included these and other data in our base¬ 
line forecasts as they are, arguably, more prone to sys¬ 
tematic uncertainties such as tracer bias and the shape 
of the non-linear power spectrum, which contains un¬ 
known baryonic feedback effects, source distribution un¬ 
certainties and multiplicative bias [II1I31I77]. However, 
a promising path would be to examine a suite of differ¬ 
ent complementary probes, each in combination with the 
CMB and BAO data, to distinguish neutrino mass from 
a non-flat or non-A model. 

For example, in Fig. we illustrate schematically how 
galaxy weak lensing can break remaining degeneracies 
between two example models. Here the CMB spectra 
and BAO distances are indistinguishable with S4-I-DESI 
data (see parameters in the Appendix). 

Their galaxy weak lensing signals, which probe the 
growth of structure during the dark energy dominated 
era, are distinct, and should be distinguishable with a 
plausible future weak lensing survey. With a 2% differ¬ 
ence in cTg, robust galaxy cluster measurements should 
also help discriminate between these models [e.g., [751175]. 


VI. DISCUSSION AND CONCLUSIONS 


We have demonstrated that an indirect detection of the 
sum of the neutrino masses should be possible with up¬ 
coming CMB and BAO surveys. In the next decade a do- 
detection should be reachable, within the ACDM-t-Em^ 
model, even in the minimal mass (Em^ = 60 meV) sce¬ 
nario. We have found that this is contingent on obtain¬ 
ing improved large-scale polarization measurements from 
the CMB, which may be the hardest experimental chal¬ 
lenge. It will also be necessary to exclude degeneracies 
with other plausible extensions such as curvature and 
dark energy. We find that allowing for these extensions 
degrades the expected neutrino mass constraint, and that 
use of other large-scale structure probes will be necessary 
to rule out other departures from ACDM. 

Our forecasts make a number of assumptions. We have 
neglected non-Gaussian terms in the CMB covariance due 
to correlation of the temperature, polarization, and lens¬ 
ing fields, which should have a small impact on immi¬ 
nent data but will get more important as noise levels 
reduce. We assume Gaussianity of the posterior distri¬ 
bution and rely on data-independence of the covariance 
matrix. These assumptions become increasingly accurate 
for the precision future experiments considered here, but 
remain an approximation. Our analysis also assumes lit¬ 
tle foreground contamination. This will be most impor¬ 
tant for the large-scale CMB signal, but is not expected 
to significantly degrade the smaller-scale lensing signal. 

We have assumed white noise in the new CMB data. 
This has not yet typically been achieved in practice from 
the ground, due to atmospheric and scan-synchronous 
emission which induces additional variance at large 
scales. However, polarization measurements are cleaner 
than temperature in this respect, and new experiments 
have sophisticated designs to modulate the polarization 
signal. Performance of current ‘Stage-2’ experiments will 
help refine noise projections for the future experiments. 
We have also assumed ideal performance of the future 
DESI BAO experiment. 

In terms of theoretical scope, beyond an evolving 
dark energy equation of state, we have not considered 
other possible new physics. For instance, theoretically- 
motivated axion contributions to the background expan¬ 
sions and perturbations might offer an alternative expla¬ 
nation for a neutrino mass detection. The cosmological 
effects are not identical [SD] , but further investigation will 
be useful. 

Finally, neutrino mass is the next beyond-ACDM pa¬ 
rameter that we know will be needed to fit data, so 
it is valuable that competitive constraints are expected 
to come from different combinations of cosmological 
datasets beyond the CMB/BAO considered here. Their 
complementarity will aid in convincingly excluding sys¬ 
tematic effects and alternative cosmological models. 
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Appendix: Forecasting parameter constraints 

We forecast posterior parameter constraints and de¬ 
generacy directions by evaluating the relevant Fisher ma¬ 
trix as defined in Eq. We assume uniform priors on all 
parameters, allowing the posterior p{6) to be replaced 
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Parameter 

Fiducial value Step size 

Q,bh^ 

0.0222 

0.0008 


0.1197 

0.0030 


1.0409 

0.0050 

10® A« 

2.196 

0.1 

Us 

0.9655 

0.010 

T 

0.078 

0.020 

(meV) 

60 

20 


0 

0.01 

Wq 

-1 

0.3 

Wa 

0 

0.6 

Neff 

3.046 

0.080 

Oe 

0.007 

0.002 


TABLE III: Fiducial values and step sizes for the numerical 
derivatives, chosen to be small enough to minimize error in the 
Taylor expansion of the two-sided derivative, while keeping 
numerical stability in the derivatives from the CAMB power 
spectra. The fiducial value of = 0.007 is at the upper 95% 
confidence level given current data m- 


likelihood into Eq. and neglecting parameter depen¬ 
dence in the power spectrum covariance matrix one ob¬ 
tains 


c«=i: 

I 


dCj^_,dCi 

de, ^ dBj ■ 


(A.3) 


From Eq. |A.1[ and applying Wick’s theorem, the covari¬ 
ance matrix for the power spectra has elements 


C{Cf,Cf) = (- 2; + 1)/ —)(gf' + <) 

+ (Cr^ + lVr^)(C'f^ + iVf^)], (A.4) 

where a, l3,^,6 G {T, E, B, Kc\ and /sky accounts for the 
loss of information due to partial sky coverage [HiiiHa. 
Noise spectra are generated for each observable given in¬ 
put noise properties such as CMB map sensitivities. We 
assume additive white-noise for the CMB: 

Nr = (AT)2 exp ( ^(^+ yFWHM \ (a.5) 

y 8 In 2 J 


with the likelihood C{6) = F{d\6,M) in Eq. The 
Fisher matrix is evaluated at fiducial parameters] the 
choice is unimportant under the assumption of a Gaus¬ 
sian posterior and data-independent covariance, although 
in practice they are chosen to match the current best-ht 
parameters of the model, as the real covariance does have 
a term that scales with the signal. Our choices are shown 
in Table UTIl 

To compute the Fisher matrix we use as observables 
the CMB power spectra, and the BAO distance measure¬ 
ments. For the CMB, we use the lensed power spectrum 
between each pair of (assumed Gaussian) fields X, Y from 
their spherical harmonic coefficients: 


cr 


1 


m—i 

'y ] 

m——t 


(A.l) 


for a G {T, E, B}, where AT {AP for polarization) is the 
map sensitivity in /iK-arcmin and 0fwhm is the beam 
width. This as an optimistic approximation: real noise- 
spectra from ground-based experiments have a dominant 
contribution from atmospheric variance at large scales 
(see e.g. Fig. 4 of Das et al. [72]). The atmosphere 
is weakly polarized, and hence the white-noise approxi¬ 
mation is better in E and B than T. The CMB lensing 
reconstruction noise is calculated using the quadratic- 
estimator formalism. As described in the main text, we 
neglect non-Gaussian terms in the power spectrum co- 
variance, and also neglect the BB spectrum as it does 
not contribute significantly to upcoming constraints and 
has a highly non-Gaussian covariance [35] . 

We add information from Baryon Acoustic Oscillation 
(BAO) experiments by computing the BAO Fisher ma¬ 
trix: 


This formula omits beam smoothing effects and the sub¬ 
traction of a noise term, which we account for below. 
The estimated power spectrum is a sum of many ran¬ 
dom variables of finite variance, and to good approxima¬ 
tion follows a Gaussian distribution. This approximation 
breaks down at large scales but does not have a signifi¬ 
cant impact on expected errors. For a full-sky survey, we 
have 

-2\nC{9) = -2'^lnp{Ce\e) 

i 

= ^ [{Ci-nercj\e){Ce-cAe))+indet{27rne)) 

I 

(A 


pBAO ^ Sr ^ 

V 


(A.6) 


where fk = f{zk) = ‘rs/dv{zk) is the sound horizon at 
photon-baryon decoupling over the volume distance 
dv to the source galaxies at redshift Zk- These real and 
forecast data are reported in Table [V| 

The total Fisher information matrix is given by the 
sum of the CMB and BAO Fisher matrices, and is in¬ 
verted to forecast parameter covariances. An alternative 
MCMC approach using simulated data can be taken to 
account for non-Gaussianity of the posterior [e.g., 1^ . 
but the Gaussian approximation is likely increasingly 
2)good as the data quality improve from Planck through 
S3 to S4. 


where Cg = (Cj^, ...) contains auto- and cross- Our forecasting code, OxFish, has been developed for 

spectra and is their covariance matrix. Inserting this this analysis and is used to forecast parameter covariance 
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Experiment /sky v/GY^t. /min /max FWHM/arcmin AT//iK-arcmin AP/^K-arcmin 

PJancJf-2015 0.44 30,44,70,100,143,217,353 2 2500 33,23,14,10,7,5,5,5 145,149,137,65,43,66,200 -,-,450,-,-,-,- 

PJanci-pol ” -,-,450,103,81,134,406 

WMAP-pol 0.74 33,41,64,94 2 1000 41,28,21,13 -,-,298,296 425,420,424,- 


TABLE IV: Specification for the Planck and WMAP experiments used in the analysis, assuming white noise properties. We 
define ‘Planck-2015’ to reproduce the constraints from Planck 2015 data; for ‘Planck-poU we use the Planck Blue Book scaling 
factors to convert to polarization. For WMAP-pol we recover an optical depth uncertainty that matches the WMAP9 data. 
When combining with S3 and S4, we include Planck-pol data across the full /sky = 0.44 at large scales {£ < across 

/sky = 0.2 in the multipole range < £ < 2500 as the useful Planck data and S3/S4 will likely not overlap completely. Our 

results are insensitive to the exact choice of this non-overlapping region size. When using WMAP-pol data we substitute it in 
at large scales {£ < /sky = 0.74. 


Experiment Redshift 

aps/dv) 

(r„/dv) 

(%) 

a{rs/dv) Ref 

6dFGRS 

0.106 

4.83 

0.0084 


SDSS MGS 

0.15 

3.87 

0.015 

m 

LOW-Z 

0.32 

2.35 

0.0023 

HU 

C-MASS 

0.57 

1.33 

0.00071 

m 

DESI 

0.15 

1.89 

0.0041 

m 


0.25 

1.26 

0.0017 



0.35 

0.98 

0.00088 



0.45 

0.80 

0.00055 



0.55 

0.68 

0.00038 



0.65 

0.60 

0.00028 



0.75 

0.52 

0.00021 



0.85 

0.51 

0.00018 



0.95 

0.56 

0.00018 



1.05 

0.59 

0.00017 



1.15 

0.60 

0.00016 



1.25 

0.57 

0.00014 



1.35 

0.66 

0.00015 



1.45 

0.75 

0.00016 



1.55 

0.95 

0.00019 



1.65 

1.48 

0.00028 



1.75 

2.28 

0.00041 



1.85 

3.03 

0.00052 



We also forecast the neutrino mass constraint from 
P15-I-BAO, finding a{Tjm^) = 103 meV. Placing the peak 
of the posterior at the fiducial = 60 meV, this cor¬ 
responds to Stoj/ < 245 meV at 95% confidence, compa¬ 
rable to the actual result of Smy < 230 meV (TH]. This 
includes JLA SNe data and an iLg prior, but these are ex¬ 
pected to have a small impact. We also agree with [2HU2H] 
on the neutrino mass constraint for the S4-I-DESI data 
combination, finding = 15 meV if we assume 

that the reionization bump is measured. 

W e u se the following parameters for the curves in 
Fig. 10 solid curve; = 0.0222, = 0.0120, 

= 3 meV, r = 0.067, IO'^Oa = 1.0417, 10®= 2.15, 
Us = 0.9647, as = 0.835, Hq = 68.0}, dashed curve; 


{flfc/i® = 0.0222, = 

T = 0.089, lO^dA = 1.0399, 
cTg = 0.822, Ho = 66.6}. 


0.0119, Emi/ 
10®A, = 2.24, 


= 117 meV, 
ris = 0.9663, 


TABLE V: Specification for current BAO-15 data (top), and 
forecast DESI data (bottom). We derive the expected frac¬ 
tional uncertainties on rsjdv for DESI from the fractional 
errors on DaIts and H{z) forecast in [22]. The absolute val¬ 
ues correspond to a ACDM model with Errii/ = 60 meV. 


matrices in one coherent python package. The code inter¬ 
faces with the CAMB code for evaluation of power spectra. 
We compare with real data or previous work where possi¬ 
ble. We construct the ‘Planck-2015’ (P15) specification, 
given in Table IV to produce constraints which match 
the ACDM uncertainties in Planck Collaboration et al. 
m, with the beam sizes and noise levels matching the 
detector sensitivities in 1551. 














